Report for Model 9

— for both jointly and separately

2024-08-15

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.06       1.17
## dh0           1.72       2.67
## dh1           1.88       2.99
## dl0           4.07       8.16
## dl1           2.50       4.81
## dl2           1.60       2.48
## dl3           2.88       5.05
## 
## Multivariate psrf
## 
## 3.64

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1403.362 23642.734
DIC3 1276.180 20946.750
PWAIC 103.750 1052.198
WAIC 1355.716 21358.294

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.02
## LDevsum       1.00       1.01
## dh0           1.20       1.59
## dh1           1.19       1.58
## dl0           1.28       1.74
## dl1           2.04       3.67
## dl2           1.16       1.42
## dl3           2.43       5.08
## 
## Multivariate psrf
## 
## 2.1

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1448.01985 23921.990
DIC3 1292.89309 21168.527
PWAIC 94.59621 1536.806
WAIC 1360.96116 21915.208